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We compute the instantaneous contributions to the spherical harmonic modes of gravitational 
waveforms from compact binary systems in general orbits up to the third post-Newtonian (PN) 
order. We further extend these results for compact binaries in quasi-elliptical orbits using the 
3PN quasi-Keplerian representation of the conserved dynamics of compact binaries in eccentric 
orbits. Using the multipolar post-Minkowskian formalism, starting from the different mass and 
current-type multipole moments, we compute the spin-weighted spherical harmonic decomposition 
of the instantaneous part of the gravitational waveform. These are terms which are functions of 
the retarded time and do not depend on the history of the binary evolution. Together with the 
hereditary part, which depends on the binary’s dynamical history, these waveforms form the basis 
for construction of accurate templates for the detection of gravitational wave signals from binaries 
moving in quasi-elliptical orbits. 

PACS numbers: 04.30.-w, 04.25.-g, 04.25.dg, 97.60.Jd 


I. INTRODUCTION 

Compact binary systems composed of neutron stars and/or black holes are one of the most promising sources for 
the second generation of earth-bound gravitational-wave (GW) detectors such as Advanced LIGO [!;) and Advanced 
Virgo jlj as well as for the proposed space-based detector eLISA Q. Detection of such systems in GW detectors relies 
on a data-analysis technique known as matched filtering which in turn requires very accurate modeling of GW signals 
from these sources dj. The compact binaries are known to have significant eccentricities when they are formed. 
However, since the GW radiation reaction effects tend to circularize the binary’s orbit Si, for most long-lived 
binary systems one can expect that their orbits would have circularized by the time they enter the sensitivity band of 
ground-based detectors. This has motivated the GW community to perform searches of GW signals from coalescing 
compact binary systems (CCBs) using circular orbit templates. 

Many astrophysical scenarios have been proposed which suggest the possible existence of close eccentric binary 
systems. 1 One such scenario may exist in the cores of dense globular clusters due to a mechanism known as the “Kozai 
mechanism” (7|. This mechanism can also come into effect in scenarios involving formation of hierarchical triples 
of supermassive black holes due to subsequent mergers of galaxies Q. Another scenario might involve formation of 
close eccentric compact binary systems in dense stellar systems like globular clusters @. Compact binaries involving 
intermediate mass black holes in globular clusters might be seen in the eLISA band with residual eccentricities of 
0T< e ^0- 2 0- Other scenarios involve formation of close eccentric compact binary systems at centers of galaxies 
0 and NS-BH binary systems which can become eccentric as a consequence of multi stage mass transfer from the 
NS to the BH 0]. In light of these possibilities it becomes necessary to compute accurate waveforms accounting for 
the eccentricity of the binary’s orbit. 

A number of investigations concerning the sensitivity of searches using circular orbit templates to detect eccentric 
binary systems have been performed in the past. The first such investigation was presented in Ref. [0 where the 
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1 Eccentric binary systems, with small orbital separations, formed through the capture process. Unlike the long-lived CCBs such systems 
are expected to enter the detector with nonzero eccentricities. 
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authors studied the loss of signal-to-noise ratio in detecting signals from binaries with residual eccentricities using 
template banks constructed with quasi-circular waveform models with leading order effects (both conservative and 
secular). They argued that even if the system has a residual eccentricity (eo 1$ 0.13 for binary system with two 1.4 Mq 
neutron stars or eo < 0.3 for a binary with two 6 Mq black holes), use of circular orbit templates will be sufficient 
to detect signals from such systems. 2 However, this result has been subsequently weakened due to two independent 
investigations BO. Both of the investigations suggest that if the eccentricity of the binary when it enters the 
sensitivity band of detector is greater than 0.1, then it will not be possible to detect such systems using circular 
orbit templates. These investigations only dealt with sources that will be seen in ground-based detectors. However, 
the capabilities of circular waveforms to detect signals from the coalescence of supermassive black holes (visible in 
the eLISA frequency band) have been investigated in [16:]. The results presented in Ref. [l|| suggest that even to 
search signals from sources with initial eccentricities of the order 10“ 4 one would need waveforms which accurately 
account for the effects of eccentricities. In addition, in a recent work, Huerta and Brown da showed that searches 
for CCBs with eccentricity > 0.05 would require eccentric template banks to avoid significant loss in the sensitivity 
of the search. Lastly, systematic errors due to the orbital eccentricity in measuring the source parameters of double 
NS systems was investigated recently by Favata [l8j which again indicated the necessity to incorporate the effects of 
eccentricity to measure the parameters of a double NS system if it has non-negligible eccentricity when detected. 


Evolution of a compact binary system can be divided into three stages: the early inspiral, late inspiral and 
merger and the final ringdown. The early inspiral phase can be very well modeled using the approximation schemes 
such as multipolar post-Minkowskian (MPM) approximation matched to post-Newtonian(PN) [19j whereas the late 
inspiral, merger and ringdown phases can be modeled using numerical relativity (NR) [20} or effective one body 
approach [2l} . In fact, it is now possible to perform numerical simulations to track the evolution of the BH binary 
systems over many inspiral orbits and the subsequent merger and ringdown phases. However, high computational 
cost of generating numerical waveforms covering the entire parameter space of coalescing binary black holes (BBHs) 
has led to the construction of hybrid waveforms (by combining PN and NR waveforms), which further are used to 
phenomenologically construct a waveform model which has sufficient overlap with the hybrid waveform [B127} . In 
addition to this, one needs to check the consistency between these two waveforms (PN & NR) in a regime where both 
of them are valid. This would not only tell us about the compatibility of the two waveforms but also would indicate 
the limits up to which PN waveforms are reliable. There have been many such investigations involving nonspinning 
and nonprecessing BBH in quasi-circular orbits f28l fiisl and quasi-eccentric orbits [39}. The need for such comparisons 
and matching of the two waveforms (PN and NR) has led to the high accuracy computations of spherical harmonic 
modes of the PN waveforms in case of CCBs moving in quasi-circular orbits |40l442| . Evidently, in order to perform 
similar comparisons for eccentric binaries, one would need high accuracy eccentric PN waveforms for such systems. 


The leading order (or Newtonian) expressions for the GW polarizations (h+ and h x ) were obtained in the context 
of spacecraft Doppler detection of GWs from an isolated compact binary in an eccentric orbit [43}. This work was 
then extended to 1PN and the next 1.5PN order in 0 - 1481 . At the 2PN order, the transverse-traceless radiation field 


(hJ T ) due to an isolated binary c omp osed of two compact stars moving in eccentric orbits was computed in [491 15 


Although the two works, [49[ and [50j, followed two different approaches, their final findings were in perfect agreement 
with each other. Under the adiabatic approximation, associated 2PN GW polarizations (/i_|_,/i x ) were obtained in 
m for the inspiral phase of binaries in quasi-eccentric orbits. Later in Refs. [52i . [53} , the method of variation of 
constants was used to compute post-adiabatic corrections (varying on the orbital time scale and 1/c 5 times smaller) 
to the secular variation due to radiation reaction. Using the 3PN generalized quasi-Keplerian representation of the 
conservative dynamics of compact binary systems with arbitrary mass ratios moving in eccentric orbits presented in 
[54j. Ref. [Ill provides the evolution of the orbital phase with relative 1PN accuracy (absolute 3.5PN). The en ergy 
and angular momentum fluxes as well as evolution of orbital elements up to 3PN order was calculated in Refs, [ool 1-571 . 
Recently, computations of the frequency domain waveforms and the orbital dynamics (both at the 2PN order) 
were presented for eccentric binaries in harmonic coordinates [58}. On the NR front, the first simulations involving 
nonspinning equal mass BBHs in bound eccentric orbits were performed in Bill and the effects of eccentricity on 
the final mass and spin were studied. Another recent work [39] presents numerical simulations for a nonspinning equal 
mass binary system with an initial eccentricity of e ~ 0.1 and compares the NR waveforms with those of the PN models. 

In this paper we present the computation of instantaneous part 3 of various modes of the waveform (h im ) for 


2 They chose a lower cutoff for the fitting factor (FF|min = 90%) corresponding to a loss in event rates of about less than 27%. 

3 The part of the gravitational radiation which depends on the state of its source at a given retarded time. 
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general orbits using the basis of spherical harmonics of spin weight -2. In addition we also specialize to the case 
of compact binaries in quasi-elliptical orbits and provide 3PN instantaneous expressions for various modes using 
3PN quasi-Keplerian representation of the conserved dynamics of compact binaries in eccentric orbits [HQ. Note 
again that investigations presented here involve only the contributions from the instantaneous terms which must be 
complemented by computations accounting for the hereditary effects. 4 Computations of hereditary parts to various 
modes of the waveform will form the basis for a companion paper [filj |. 

This paper is organized in the following manner. In Sec.[TT]we first introduce general formulas for spherical harmonic 
modes of the gravitational waveform, h fm , in terms of the radiative mass and current multipole moments, U im and 
V tm . Section Till Al recalls some of the important aspects of the MPM-PN formalism and lists various inputs that 
are needed for computing 3PN expressions for various modes. These inputs involve relations connecting radiative 
moments to source moments, expressions for various source multipole moments for an isolated compact binary system 
and equations of motion. In Sec. II VI we provide our results related to the instantaneous part of the spherical harmonic 
modes of the waveform for a nonspinning compact binary system in terms of variables that describe the radiation from 
a generic compact binary. We find that these expressions are quite large and run over several pages. Keeping this in 
mind we choose to list only the dominant mode (/ 122 ) in the main text of the paper and provide all the relevant modes 
contributing to the 3PN waveform in a separate file readable in MATHEMATICA (Hlm-GenOrb.m) that will be made 
available on the journal web page as Supplemental Material [62] along with the paper. In Sec. [V] we specialize to the 
case of CCBs moving in quasi-elliptical orbits and provide the corresponding expression for the dominant mode, h 22 , 
in terms of the time-eccentricity et, a PN parameter related to the orbital frequency x and the eccentric anomaly u. 
Similar to the general orbit case, in the case of CCBs in quasi-elliptical orbits, expressions for all the relevant modes 
contributing to the 3PN waveform will be listed in a separate file (Hlm-EllOrb.m). Finally in Sec. IVII we conclude 
the paper by providing a brief summary of our results and the future plans. 


II. SPHERICAL HARMONIC MODES OF THE GRAVITATIONAL WAVEFORM 


For an isolated source of GWs, the spherical harmonic modes of the waveform (h em ), in terms of the radiative 
mass-type (U em ) and current-type multipole moments ( V tm ) 1(33|. are given as 


h lm = - 


G 


V2Rc e + 2 


U em _ ly 

c 


lm 


( 2 . 1 ) 


where R is the distance of the source in radiative coordinates, G is Newton’s gravitational constant and c is the speed 
of the light. The radiative multipole moments, U im and V fm , appearing above are related to the symmetric trace-free 
(STF) radiative moments Ul and Vl as 
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Here denote STF tensors which connect the usual basis of spherical harmonics Y (m (€),<!>) to the set of STF 
tensors N (L ) = N( h ■ ■ ■ N it ) as 5 
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4 The part of the gravitational radiation which depends on the entire dynamical history of the source and is complementary to the 
instantaneous part of the radiation. 

5 Here L = represents a multi-index composed of l spatial indices and the angular brackets (()) surrounding indices denote 

symmetric trace-free projections. 
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In the above, N = X./R is a unit vector pointing towards the detector along the line joining the source to the 
detector. For instance, if the binary’s plane is the x-y plane then N, in terms of angles (0, $) giving the location of 
the binary, can be given as 

N = sin 0 cos $ e x + sin 0 sin $ e y + cos Qe z . (2-4) 

The STF tensorial coefficients o^ m in terms of ■ ■ ■ N if \ and Y lm (Q 1 4>) can be written as 6 

a{ m = JdnN (L) Y em , (2.5) 

where the usual basis of spherical harmonics is given as 
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( 2 . 6 ) 


It is important to note that for nonspinning binaries, there exists a mode separation as pointed out in Ref. ]4l| and 
explicitly shown in Ref. [42]. The mode h lm is completely determined by mass-type radiative multipole moment 
(U im ) when l + m is even, and by current-type radiative multipole moment (V^ m ) when t + m is odd. This allows us 
to write for various modes 


h tm = - G utm if e + m is even, (2.7a) 

V2Rc*+ 2 K ’ 

h* m = — - V £m if £ + m is odd. (2.7b) 

V2Rc e + 3 V ' 


III. INPUTS FOR COMPUTING THE 3PN WAVEFORM 
A. Relations connecting the radiative moments to the source moments 

In the MPM-PN formalism [H, Ifui - liki ] , the radiative multipole moments (Ul, Vl) are first written in terms of 
two sets of canonical moments (-Ml, Sl), which in turn are expressed in terms of six sets of source moments (/l, 
Jl , Wl, Yl , Xl, Zl). Relations connecting radiative moments to the canonical moments and those connecting the 
canonical moments to the source moments, with the PN accuracy desired for the waveform computations at the 3PN 
order, have been established and have been listed in Ref. [40] (see Eqs. (5.4)-(5.8) and Eqs. (5.9)-(5.11) there). Using 
these inputs we can parametrize the set of radiative moments (and, hence, the modes) in terms of source multipole 
moments. Below we list all the relevant radiative multipole moments in terms of the source multipole moments with 
PN accuracy desired for the present work. Further, these expressions can be decomposed into two parts namely the 
instantaneous contribution and the hereditary contribution. 

The only radiative moment required at the 3PN order is the one related to the mass quadrupole (Uij) and is given 
by 


tt _ rrinst tt hered 

u ij — u ij ' u ij j 


(3.1) 


where the instantaneous and hereditary parts in terms of the source multipole moments read 
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In the above, the quantity M represents the mass monopole moment or the Arnowitt, Deser and Misner (ADM) 
mass of the source. The constant tq appearing in the above integrals is related to an arbitrary length scale tq by 
r 0 = To/c and was originally introduced in the MPM formalism. Note that, numbers in the parenthesis (appearing 
as superscripts of the source moments) denote the p th time derivatives. The Levi-Civita tensor is denoted by £ijk, 
such that £123 = +1 and 0(l/c 7 ) indicates that we ignore contributions of order 3.5PN and higher. 


As may be seen from the above, computing the instantaneous part requires source multipole moments given at a 
retarded time U. On the other hand, the hereditary part involves integrals over time and would require the knowledge 
of the source multipolc moments at any instant of time before U in the past dynamical history of the source. Further, 
the hereditary terms are of two kinds: those with and without the logarithmic factors (see Eq. (I3.2bl) above). The 
first integral appearing in Eq. (I3.2bl) (the one with the logarithmic kernel inside) is called the “tail-integral” and the 
one in the last line is called the “tail-of-tail” integral whereas the integral without the logarithmic factor (in the 
first line) is known as the “memory” integral. This paper only focuses on computing the instantaneous contribution 
to various modes of gravitational waveforms and the computation of hereditary contributions shall be discussed 
elsewhere 111 - 

Moments required with 2.5PN accuracy are the mass octupole moment Uijk and the current quadrupole moment 
Vij. The mass octupole moment U^k is given as 
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The current quadrupole moment VA is given as 
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where V) lnst and V)') ered in terms of source multipole moments read 
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At the 2PN order the required moments are Uijki and Vijk- The moment Uijki is given by 
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where 14™ st and V^ red are given in terms of the source multipole moments as 
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The moments required at the 1.5PN order are Uijkim and Viju- The mass-type moment Uijkim is given as 
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The current-type moment Vy.( is given by 

where l/yf and V^|j ed in terms of the source multipole moments read 
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Other mass-type moments f/^ contributing to 3PN waveform are given as 

U L = Uft + t/£ ered , 

where Z7™ st and t/ dered are related to the source multipole moments as 

U i ft t (U)=ljp(U) + 0 
Uft ed {U) = O 

Other current-type moments Vl contributing to 3PN waveform are given as 
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where finally V)! nst and Vft ed in terms of source multipole moments read 
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B. Source multipole moments in general dynamical variables 

What we need next are expressions for various source multipole moments with the PN accuracy sufficient for 
the present computation. Expressions for various multipole moments presented here are generalizations of related 
circular orbit expressions presented in [ioj, [69:j to the case of general orbits and have been computed using the 
methods presented in [651, @1 ■ We skip all the details of the computation and list the final expressions for the source 
multipole moments related to a source composed of two nonspinning compact objects moving in general orbits. 

The only moment required here with 3PN accuracy is the mass quadrupole, Iij , which for CCBs in general orbit 
was computed in Ref. [ 66 | and listed in Ref. 0] in standard harmonic (SH) coordinates . 7 As was argued in Ref. [bfj , 
though the use of SH coordinate is useful in performing algebraic checks on PN computations, quantities when 
expressed in these coordinates involve some gauge-dependent logarithmic terms and are not suitable for numerical 
calculations. It was suggested in Ref. [56} that such logarithms can be transformed away by using some coordinate 
transformations. They showed how the use of a modified harmonic (MH) coordinate system (or alternatively an ADM 
coordinate system) removes these logarithms. We skip the details related to those transformations and directly write 
the expression for the mass quadrupole moment in MH coordinates. In MH coordinates, to 3PN accuracy, hj reads 
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7 Note that Ref. [56| lists explicit expressions for all the source multipole moments for binaries in general orbits needed for computing 
3PN energy flux. 
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(3.20a) 


(3.20b) 


(3.20c) 


In the above, x t and Vi denote the binary’s relative separation and relative velocity, respectively, whereas r denotes the 
radial velocity. As we see, the above expression still has a dependence on some logarithms (log[ro]), where the quantity 
ro is related to the arbitrary constant tq appearing in tail integrals by tq = ro /c. It has been argued and shown that it 
disappears from all the physical quantities like the radiation field at infinity and the far-zone energy flux [40ll4ll . 5£ 


The expression for mass octupole, Iijk, at 2.5PN order reads 


hjk = — ism A 

r 2 
6 


„ 56 is G 2 m 2 . 
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(3.21) 
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where, 
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G 2 to 2 / 47 1591i/ 235z/ 2 \ 4 /257 7319 z/ 
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The remaining mass-type source multipole moments with PN accuracy required in the present work are 
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(3.22a) 

(3.22b) 

(3.22c) 

(3.22d) 


(3.23a) 

(3.23b) 

(3.23c) 

(3.23d) 

(3.23e) 
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The current quadrupole moment is needed at 2.5PN order and given as 


where, 
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Other current-type source multipolc moments with PN accuracies sufficient for present calculations read 
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Jijklmno ^ m £ab(i‘Ejklmno) •Ea'Vb (l "7 V 14 IS 7 IS (D . 


(3.26e) 

The required gauge moments, the monopolar moment W, and two dipolar moments Wi and P) are finally given by 

(3.27a) 
(3.27b) 


W=-ismrf + o(— ) , 


1 


Wi = ism A r 2 


v l - 3 - x l 
r 


■o\ 4 


Yj = — ism A 
5 


lGm 


1 


x l — v 2 x l — - rr v * 
2 r 2 2 


+ g 


(3.27c) 


C. The post-Newtonian compact binary dynamics 


Since relations connecting the radiative multipole moment to the source multipolc moment involve time derivatives 
of the source multipole moments, computations of various modes will require a knowledge of the equations of motion 
(EOM) with the PN accuracy with which one wants to compute various modes. Before we write expressions for 
the EOM, with the PN accuracy required for the present work, let us recall the definitions of various dynamical 
variables as well as some other related results which will be used in calculations performed here. The binary’s relative 
separation, x 1 is given by 

% l = y\ - V l 2 = r n% (3.28) 

with r = |x|, where x is the relative separation vector. Here, y\ and y 2 are position vectors of the individual 
components of the binary and n 1 are components of the unit vector, n, along the relative separation vector. 

For the relative velocity and relative acceleration we have 


v l = - and r = n ■ v, 

dt 

_ dv l _ d 2 x i 
a * ~ ~dt ~ W 

In addition we also would need expressions for v and r, which can be given as 


(3.29) 

(3.30) 


a • v 


v 


r = — [(i> 2 — r 2 ) + a • x] . 


(3.31) 

(3.32) 


with v = |v|. 

While computing the time derivatives of the source multipole moments, whenever quantities like a 1 or v or f appear, 
they are consistently replaced by their expressions in terms of variables related to the position and velocities (r, r, v). 
Computations of various modes at the 3PN order would require the knowledge of the 3PN EOM governing the 
compact binary dynamics. EOM associated with SH coordinates at 3PN order for a system of two compact objects 
moving in general orbits is available in the literature [U, |63, [7l| - l73| and are given in terms of the variables related to 
the position and the velocity of individual constituents of the binary. Since we are using expressions for the source 
multipole moments in the center-of-mass frame of the system, we need the EOM reduced to the CM frame. The 
3PN accurate expression for relative acceleration, reduced to the CM frame, in SH coordinates, were obtained in M- 
However, as discussed in the previous section (about using MH or ADM coordinates instead of SH coordinates), we 
wish to use EOM reduced to CM frame associated with MH coordinate which is given in [75| and takes the following 
form, 
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(3.33) 
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where, 
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(3.34a) 


(3.34b) 


We now have all the inputs which are needed to compute the instantaneous expressions for various spherical 
harmonic modes (h em ) associated with 3PN gravitational waveforms of GW signals from CCBs moving in general 
orbits. With this motivation, we shall proceed towards the next section where we shall present our results. 


IV. INSTANTANEOUS TERMS IN THE 3PN GRAVITATIONAL WAVEFORM FOR CCBS IN 

GENERAL ORBITS 

Combing Eq. (12.21) and Eq. (12.71) we can write the instantaneous part of various modes as 


h 
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G 4- I (£+ l)(l + 2) lm - t 

V2r c ^£\ V m-v 

if i + m is even, 
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V2Rc e+3 £'■ y 2(£ +!)(£— !) L L 


\f£ + m is odd. 


(4.1) 


(4.2) 


Relations connecting the instantaneous part of STF radiative moments {U™ Bt and V™ st ) to the source multipole 
moments have been listed in the previous section. This allows one to write the instantaneous part of various modes 
te) in terms of the source multipole moments. With expressions for the source multipole moments for CCBs 
moving in general orbits and their relevant time derivatives, one can write expressions for various modes in terms of 
dynamical variables related to the position and velocity (r, r, <j>, v). 8 


8 Alternatively, one can also compute various modes associated with the gravitational waveform using polarization waveforms (see Sec. 
II and IX of Ref. [40ll for the details). 
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Again, since v 2 = r 2 + r 2 (jr , we can write various inodes of the waveform in terms of the dynamical variables, 
namely, the radial separation (r), radial velocity (r), orbital phase and the angular velocity (</>) . The structure 
of h em reads 
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(4.3) 


For the dominant mode {l = 2, to = 2), with 3PN accuracy, various PN pieces of the coefficient H 22 st read 
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(4.4c) 


(4.4d) 
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(4.4e) 


As mentioned in Sec. |U since expressions for various modes are very large and would run over several pages, we have 
decided to provide an additional file (Hlm-GenOrb.m) containing expressions for all the modes which will be made 
available along with the paper. Finally, circular-orbit limit of the instantaneous hg m can be obtained by replacing 
related expressions for cj){= w), f and r given in Sec. IV of [40j. 

Note that he m can directly be used to write the polarization waveforms (/z+, h x ) using the standard decomposition 
of h + and h x in terms of spherical harmonic modes of spin weight -2 given in Ref. [40L |41|| . 


+ 00 t 

h + ^ih x = J2 E 

t=2 m——i 


(4.5) 


where Y_f™’s (the spin-weighted spherical harmonics of weight —2) are functions of the spherical angles (0, <F) defining 
the binary’s location and given as 


y tm _ 

-2 — 




47T 


4 m (0)e* 


k2 


d{ m = E 


{-) k y/{£ + m)!(£ - m)\{£ + 2)!(£ - 2)! ( 0\ 


k—k i 


k\ (k — m + 2)!(^ + m — k)\(i — k — 2)! 


cos ■ 


■ 


2l-\-m—2k—2 / q \ 2k— m+2 

sin ■ 


(4.6a) 

(4.6b) 


with k\ = max(0, m — 2) and &2 = min(£ + m ,i — 2). 


Further, the polarization waveform can be used to write the transverse-traceless part of the radiation field ( hij ) by 
using the following relation [i(i|, EH 


= 2 ( h + e tj + h x e ij) > 


(4.7) 
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where 


4 = \ (■ PiPj - QiQj) , (4.8a) 

4 = l(PiQi + PjQi) ■ (4-8b) 

Here P and Q are the two unit polarization vectors and they have been chosen following the convention used in [4(| . 


V. QUASI-KEPLERIAN REPRESENTATION 

In the previous section, we presented general orbit expressions for the dominant mode of the gravitational waveform 
from a compact binary system. In this section we aim to specialize to the case of compact binary systems in quasi¬ 
elliptical orbits. Expressions for various h(, m describing the radiation from binaries in quasi-elliptical orbits can simply 
be obtained by using relations connecting generic dynamical variable r, r, <j> and <j> to a set of parameters associated 
with elliptical orbits in general orbit expressions for various modes. Such relations can be established by using the 
generalized quasi-Keplerian (QK) representation of the conservative dynamics of the binary moving in eccentric or¬ 
bits, which indeed is available to us due to the work of Memmesheimer, Gopakumar and Schafer (hereafter MGS) [54] ]. 

The QK representation was first introduced by Damour and Deruelle [76| and dealt with the binary dynamics 
at 1PN order. The generalized QK representation at 2PN order in ADM-type coordinates was given in Refs. 0- 
|79| . MGS provides the 3PN generalized QK representation in both the ADM and MH coordinates, which involve 
expressions of the orbital elements associated with the orbit of the binary in terms of the conserved energy and orbital 
angular momentum of the binary. Before we get into the details of the parametrization, we first summarize equations 
describing the radial and angular motion of the binary in terms of various orbital elements associated with elliptical 
orbits (see Refs. [53], [56| for details). In the parametric form, the radial separation, r, is given by 

r = a r (1 — e r cos u) , (5-1) 

where a r is the semi-major axis of the orbit and e r is the eccentricity of the orbit (both labeled after the radial 
coordinate, r). The quantity u is called eccentric anomaly and at the 3PN order it is related to the mean anomaly 
( l ) by the relation 


l = u — e t sin u + f t sin V + gt [V — u] + i t sin 2V + h t sin 3V. 

The orbital phase, cj) , at the 3PN order reads 

cj) = 4>p + K [V + frj, sin2K + g^ sin3K + sin4K 
+h ( p sin5K], 


(5.2) 


(5.3) 


where tpp is the initial phase at the first passage of the periastron and V is the true anomaly that takes the form 

( 1 ) 


V =V(u, e^) = 2 arctan 


1 + 6cf> 

1 - erf, 


tan I 


(5.4) 


Also, the mean anomaly, l, is related to the time as 


l = n(t — tp) 


(5.5) 


where tp is the instant of the first passage at the periastron and n = 2-k/P is the mean motion with P being the 
orbital period. 

In addition to this, expressions for the radial and angular velocity can be given as 


/ dlY 1 dl 


r — a r e r sin u —- — 


\du J 1 9f ’ 

cj) = K (1 + 2 cos 2V + 3 g^ cos 3K + 4 cos 4V 


(5.6a) 


c , . dV f dl\ -1 dl 

+5/VCOS510— (-) 


(5.6b) 
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It may be seen from Eqs. (15.21) . (15.41) and (15.51) that 


dl 


dt 

dl_ 

du 


= n, 


= 1 - e t cos u + ft 


, r dV .dV 


ay _ (i-e/) 1/2 
(1 — cosu) 


(5.7a) 


(5.7b) 


(5.7c) 


In the above, and et denote eccentricities related to the coordinates <j> and t, respectively. K is related to the 
advance of the periastron per orbit and is given by K = $/(27r), where <f> is the angle of return to the periastron. In 
this parametrization, f t , g t and g$ contribute both at 2PN and 3PN order whereas it, i^, ht and contribute 
only at the 3PN order (see Ref. [56] for related details). 


Once we have written equations connecting the generic dynamical variables (r, r, <j>, and <j>) to the orbital elements 
of the elliptical orbit, we can use inputs from MGS to express them in terms of a suitable set of parameters of 
our choice. The main result of MGS is that, it provides 3PN accurate expressions for various orbital elements 
(a r ,et,e r ,e q j...) associated with the elliptical orbits in terms of the corresponding conserved energy per unit reduced 
mass ( E ) and the parameter h , related to the reduced angular momentum (J), by h = J/Gm . 9 Using these relations 
one can express the dynamical variables r, r, <f> and <j> in terms of E, h and u. Here one should note that this 
is not the only way in which the orbital dynamics can be parametrized. In fact, one can re-express E and h in 
terms of any of the two orbital elements to write equations describing the orbital motion of the binary; however 
a parametrization involving gauge invariant parameters is sometimes preferred as such parametrization is suitable 
for making comparisons with related numerical results. This led MGS to use a parametrization involving n and 
K = $/(2 7r) (both are independent of the coordinate system used when expressed in terms of E and h in order to 
describe the orbital motion of the binary in elliptical orbits). 10 Here, n is the mean motion and K = $>/(2n) denotes 
the angle of the advance of the periastron per orbital revolution. In a work related to the phasing of the GWs 
from inspiralling compact binary in elliptical orbit due to Damour, Gopakumar and Iyer [52| . the orbital dynamics 
has been described using n and et as parametrizing variables. Following the conventions of j52j . Konigsdofer and 
Gopakumar [53| provided the 3PN accurate expressions for r, r, <fi and (j) in terms of n , et and u. 11 Reference [5? > 
makes an alternative choice of parametrization in terms of variables x and et, where x is related to the orbital 
frequency to by, x = (Gmoj/c 3 ) 2 ^ 3 , and is independent of the choice of the coordinate system used. 12 The choice 
{x,et} as parametrizing variables leads to expressions which can be reduced to those related to the circular orbit 
case (et —> 0) which uses x as the expansion parameter. In addition to this, in another related work Hinder et al. [39] 
compared the two parameterizations and concluded that the choice of x as compared to n provides better agreement 
with NR results. They separately discuss these two PN models (based on the choice of parametrization as x or n) 
and call them x-model and n-model. Note that although [39] uses 3PN QK representation describing the conserved 
dynamics it used only 2PN accurate dissipative dynamics and, hence, when used to obtain quasi-circular orbit limit 
of the orbital phase it shows some deviations from related standard results as was pointed out in Ref. 
discussed below our results must be coupled with 3PN evolution equations for QK variables given in Ref. 


hence, would be more suitable for obtaining quasi-circular limits of the amplitude and phase of various modes. 


8C|. As 
57j and, 


Following the arguments presented above we choose to parametrize the dynamical variables {r, r, 4>, (/>} in terms 
of the QK variables {x,et,u} which can further be used to write expressions for various h( m in terms of the QK 
variables. Note that Ref. [39] already lists the 3PN expressions for r and </> in terms of {x,et,u}', however, they 
numerically differentiate r to obtain r and perform numerical integration of (j> to obtain (f> in order to avoid the use of 
long and complicated expressions for r and <j> in their numerical code. As mentioned above, Ref. [53| lists 3PN accurate 
expressions for {r, r, (j), <j)} in terms of the QK variables {n,et,u} in MH coordinates (see Eqs. (23)-(27) there). In 


9 In Ref. |56| . which uses results obtained in MGS, orbital elements are expressed in terms of the parameters {e,j} instead of { E. h \, 
where e and j are defined as e = —2 E / c 2 and j = —2 E h 2 . 

10 In fact, MGS uses xmgs = (G mn/c 3 ) 2 ^ 3 and the parameter k' = (K — l)/3. 

11 In fact, [53 | is the extension of (53 | and discusses the 3PN conservative dynamics of the binaries in elliptical orbits. 

12 Parameters n and ui are related by n = Kuj, and K has been defined above. 
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fact, they use a variable £ related to n by £ = {Gmn/ c 3 ). Hence to obtain related expressions parametrized in terms 
of {x, e t ,u}, all we need to know is how £ is related to x and e t . The relation between £ and our QK variables x and 
et with 3PN accuracy in MH coordinates is given as 


£ = 


+ x' 


r 3/2 


(l-e t 2 ) 3 
3 

2 + " 


1 — 3e 2 + 3ej) — e® + x (—3 + 6e 2 — 3ej)) + ; 


9 ( 33 v 

-2 +7i/+ (~t - 2 


457 123tt 2 

"1 32 


— 7u 2 


267 


279 123tU 


128 


— 40z/ 


e t 


51 13 v 

T ~~ — 

39 55i/ 65^ 2 

T + ~1 8~ 


e t 


+\]l-e 2 (-15 + 61/ + (-30 + 12^)e 2 ) 


(5.8) 


Using the above in Eqs. (23)-(27) of Ref. [53} one can write expressions for r, r, (f> and (j> in terms of the variables 
{ x,et,u } and then use them to obtain the expressions for the spherical harmonic modes of the waveform for 
quasi-elliptical orbits. 


We now are in a position to use the expressions for r, r, <f> and <f> in terms of the variables { x , et, it} to reexpress the 
instantaneous part of spin-weighted spherical harmonic modes in terms of parameters x, et and it. Just like 

the expressions for various h^ m in general orbits we find that the expressions for various modes in QK representation 
is too large to be given in the main text of the paper. Hence we will just provide the expression for the dominant 
mode here and list all the relevant modes in an additional file (Hlm-EllOrb.m). In addition, since even the /122 
expression runs over many pages we provide the PN structure of /122 in the main text of the paper and list ex¬ 
plicit expressions for various PN pieces in Appendix A to maintain the flow of discussion in the main text of the paper. 

The structure of various modes, remains the same as in Eg. (14.31) : however, now the coefficients H?™ t and q i 
are functions of the parameters {x, et,u}, 


uZm _ AG mUX pK fjf.m 

inst “ c 2 R V 5 inst ' 


(5.9) 


The instantaneous part of the dominant mode reads 


h 


22 

inst 


4 Gmux 
c 2 R 



(5.10) 


with 


T T - 

t/i 


22 


— TT 22 
ii Newt 


TT 22 
^lPN 


TT 22 
^2PN 


TT 22 
il 2.5PN 


TT 22 
^3PN 


(5.11) 


where various PN pieces appearing in Eq. (15.111) are listed in Appendix A. Before we move to the concluding section, 
we would like to make a few remarks about the results presented here. 


We observe logarithmic dependences on the arbitrary length scale tq in Eq. (jAlol) through the quantity xo which 
is related to ip by xq = (Gm/ c 2 tq). This dependence is due to the presence of terms involving logarithms of tq in 
the expression for the mass quadrupole moment (Rj) given by Eqs. (I3.19D - (I3.20I) . As was discussed in Sec. IIIIBI 
we would expect such dependences to disappear from final expression for various modes (for instance see Ref. [40} 
which lists h( m for binaries in quasi-circular orbits). As has been observed in Refs. |40l. PTII. MM , it turns out the 
hereditary contribution has equal and opposite dependences on the arbitrary length scale r$ and cancels out from 
the final expression. Since we do not provide hereditary contributions in this paper, such cancellation can not be 
shown here. However, we explicitly show this cancellation in [6l| which deals with computation of hereditary effects 
in various modes at 3PN order. 


In this paper we used the 3PN accurate QK representation which describes the conserved dynamics of CCBs in 
eccentric orbits to obtain various modes in terms of QK variables (x, et). However, it should be noted that the 
parameters x, et evolve with time over radiation reaction time scales and these secular effects start to show at 2.5PN 
order 0, |57}. Hence, in order to correctly account for the reactive dynamics of the binary at 3PN order our results 
should always be coupled with equations describing secular time-evolution of x and et . Equations describing the 
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secular evolution of orbital elements with relative 3PN accuracies were presented in Ref. [57} (see Sec. VI there). 
Evolution equations (due to instantaneous terms in energy and angular momentum loss) for orbital frequency (dw/dt 
(related to x by u> = (c 3 x 3 ^ 2 /Gm)) and for time-eccentricity parameter ( det/dt ) in terms of x and e* have been listed 
as Eqs. (6.14)-(6.15) and Eqs. (6.18)-(6.19), respectively, in Ref. [57} . 


VI. SUMMARY AND CONCLUDING REMARKS 

In this paper we presented computations of the instantaneous contributions to all the relevant modes of the 3PN 
accurate gravitational waveform of the GW signal from nonspinning coalescing compact binaries in general orbits. 
The expression for the instantaneous part of the dominant mode in terms of the variables r, r, (f> and <j), has 

been given by Eqs. (14.3D - (14.4I) above, whereas expressions for other subdominant modes (along with the dominant 
mode) have been listed in a separate file (Hlm-GenOrb.m) that is being made available along with the paper. 
Next, we specialized to the case of CCBs in quasi-elliptical orbits using the 3PN quasi-Keplerian representation 
of the conserved dynamics of compact binaries in eccentric orbits in Sec. El Related 3PN accurate expressions for 
the instantaneous part of the dominant mode in terms of the variables, namely, the time-eccentricity et, a 

PN parameter x and eccentric anomaly u, is given by Eqs. (15.10l) - (15.11l) and Eq. (1A1I) . The expressions for other 
sub-dominant modes (along with the dominant mode) have been listed in a separate file (Hlm-EllOrb.m) that is 
being made available along with the paper. 

We once again remind the readers that the results presented here only account for the contributions from the 
instantaneous terms in the waveform which must be complemented by computations accounting for the hereditary 
effects. Our investigations suggest that it is not possible to provide closed-form analytical expressions for the 
hereditary terms for binaries moving in general orbits. Moreover, even for the special case of CCBs in quasi-elliptical 
orbits it may not be possible to have closed-form analytical expressions for the hereditary terms valid for systems 
with arbitrary eccentricities. However, we find that such computations can be performed assuming an expansion 
in the eccentricity parameter (e*) (6f}. Unlike the results presented in this paper which can be applied to a binary 
with arbitrary eccentricity, results of can only be applied to systems with small eccentricities. However, the 

positive side of the work is that we shall have complete 3PN analytical expression for the waveform for binaries in 
quasi-elliptical orbits that can be used for comparison with related numerical relativity results, which is one of the 
motivations for high PN order computations of the gravitational waveforms. In addition, with complete waveforms 
at hand one would be able to write complete polarization waveforms (as discussed in Sec. m which would be useful 
for data-analysis purposes. 
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Appendix A: Various PN pieces of the coefficient H 22 


associated with the dominant mode h 22 
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+45824486400tt 2 ) - 7296998400t/ 2 - 18373747200t/ 3 ) + e 6 t (-10724639232 


+5206671360 log (1 - e t cos it) - 5206671360 log — + v (-375209879040 


V^o/ 

-4909766400tt 2 ) + 288255360000t/ 2 + 3982937600t/ 3 ) + ef (12237465600 
+24072192000i/+ 7993420800t/ 2 + 1070822400t/ 3 )) + (e t cos it) 6 (-411239646720 
+v (2829529175040- 74464790400tt 2 ) - 207633139200t/ 2 + 3495795200t/ 3 
+e? (-1455591651840+ v (2185925806080 - 4091472000tt 2 ) - 724992192000t/ 2 
-10543577600t/ 3 ) + ej (9839024640 + v (-146180229120 - 7364649600tt 2 ) 
+15925593600t/ 2 + 10655961600t/ 3 ) + e\ (-1425415680 + 162493286400t/ 
-69507648000t/ 2 - 3664371200t/ 3 ) + e s t (-1385856000 - 6721920000t/ 
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+422092800i/ 2 + 56192000i/ 3 )) + (e t cos u) 7 (132690700800 + v (-909674841600 
+22912243200tt 2 ) + 71153510400;' 2 - 196672000;/ 3 + e\ (518098291200 
+v (-814416422400 + 5728060800tt 2 ) + 246203596800;' 2 + 590016000;/ 3 ) 

+ej (-26003980800+ 72602611200i'+ 9760665600^ - 590016000;/ 3 ) 

+e t 6 (-4850496000- 23526720000^+ 1477324800^ + 196672000^)) 

+ (e* cos u ) 8 (—18955814400 + v (129953548800 - 3273177600tt 2 ) - 10164787200i' 2 
+28096000i/ 3 + e\ (-74014041600 + v (116345203200 - 818294400tt 2 ) 
-35171942400i/ 2 - 84288000^) + e\ (3714854400- 10371801600;' 

-1394380800i/ 2 + 84288000i/ 3 ) + (692928000 + 3360960000;' - 211046400;' 2 

—28096000;/ 3 )) + i (e t sin u) (74511360000 + v (-225308160000+ 2727648000tt 2 ) 
+13199155200;' 2 + e 2 t (-387763200000+ v (693492940800 - 7091884800tt 2 ) 

—37407744000;' 2 ) + ej (417871872000 + v (-487101542400 + 4364236800tt 2 ) 
+729907200;' 2 ) + e\ (-104620032000 + 18916761600^ + 23478681600z/ 2 ) 

+ (e t cosw) (—241173504000 + v (1018828800000 - 13092710400tt 2 ) - 56932761600;' 2 
+e? (1711176192000+ v (-3424055500800 + 36004953600tt 2 ) + 194885222400;' 2 ) 

+ef (—1894109184000 + v (2447865446400- 22912243200?r 2 ) - 12043468800;' 2 ) 

+e? (424106496000 — 42638745600^ — 125908992000;' 2 )) + (e t cos u ) 2 (196314624000 
+v (-1765980057600+ 24548832000tt 2 ) + 89961062400;' 2 + e 2 t (-3063785472000 
+v (6872258764800- 73646496000tt 2 ) - 423163699200;' 2 ) + (3526820352000 

+v (-5087179468800 + 49097664000tt 2 ) + 59122483200;' 2 ) + e? (-659349504000 
— 19099238400^ + 274080153600;' 2 )) + (e t cosw) 3 (136249344000 + v (1381795430400 
-21821184000tt 2 ) - 54499737600;' 2 + e 2 t (2898643968000 + v (-7118338867200 
+76374144000tt 2 ) + 491957452800;' 2 ) + e\ (-3523627008000 + v (5612094259200 
-54552960000tt 2 ) - 130653388800;' 2 ) +e Q t (488733696000+ 124449177600^ 

—306804326400;' 2 )) + (e t cos u) 4 (-292571136000 + v (-358171545600 
+8182944000tt 2 ) - 4926873600;' 2 + e\ (-1608532992000+ v (3959624908800 
-40914720000tt 2 ) - 324626227200^) + e 4 t (2082972672000 + v (-3485337292800 
+32731776000tt 2 ) + 146711347200;' 2 ) + e\ (-181868544000- 116116070400^ 
+182841753600;' 2 )) + (e t cos u) 5 (149935104000 - 107418009600^ + 18977587200;' 2 
+e 2 t (546974208000 + v (-1093948416000+ 9819532800tt 2 ) + 116055244800;' 2 ) 

+e 4 (—739031040000+ v (1166391705600- 9819532800tt 2 ) - 82114560000;' 2 ) 

+e? (42121728000 + 34974720000^ - 52918272000;' 2 )) + (e t cos u ) 6 (-23265792000 
+v (56253542400 - 545529600tt 2 ) - 5778432000;' 2 + e 2 t (-96712704000 
+v (110966169600- 545529600tt 2 ) - 17700249600^) + e? (129102336000 
+v (-166733107200+ 1091059200tt 2 ) + 18247680000;' 2 ) + e\ (-9123840000 

—486604800^ + 5231001600;' 2 ))) + \ j \- e \ (-39536640000 + v (-37022515200 
+1091059200tt 2 ) - 4866048000;/ 2 + e \ (-183389184000 + v (472087756800 
-5455296000tt 2 ) - 25911705600;' 2 ) + (332107776000 + v (-469492531200 

+4364236800tt 2 ) + 12773376000;' 2 ) + e \ (-109181952000 + 34427289600^ 
+18004377600;' 2 ) + (e t cosm) (285956352000 + v (94893004800 - 4637001600tt 2 ) 
+26793676800;' 2 + e 2 t (809740800000 + v (-2359414886400+ 27549244800tt 2 ) 
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+142483968000;/ 2 ) + ef (-1500187392000 + v (2355425740800 - 22912243200tt 2 ) 

—65813299200;/ 2 ) + e\ (442810368000- 106231910400;/- 103464345600;/ 2 )) 

+ (e t cos u) 2 (-702383616000 +1 / (-158623027200 + 8728473600tt 2 ) - 32176742400^ 
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+e e t (-688545792000 + 80168140800;/ + 239044608000;/ 2 )) + (e t cos u ) 3 (714016512000 
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+u (5509679155200- 54552960000tt 2 ) - 189593395200;/ 2 ) + e\ (509718528000 
+53100748800;/ — 281622528000;/ 2 )) + (e t cos u ) 4 (-163012608000 
+i/ (-1666905292800 + 19093536000 tt 2 ) + 237828096000;/ 2 + e 2 (-3119440896000 
+v (5801890406400- 51825312000tt 2 ) - 524012544000;/ 2 ) + ej (2129504256000 
+i/(-3504101990400+32731776000 tt 2 ) + 111006720000^) + e\ (-188255232000 
—94401331200;/ + 175177728000;/ 2 )) + (e t cos u ) 5 (-234862848000 
+v (2089780070400- 21548419200tt 2 ) - 284815872000;/ 2 + e\ (2485790208000 
+v (-3852277862400+ 31367952000tt 2 ) + 325021593600;/ 2 ) + e\ (-952300800000 
+v (1192592332800 - 98195328007T 2 ) + 12165120000;/ 2 ) + e\ (42577920000 
+33423667200;/ - 52370841600;/ 2 )) + (e t cos u ) 6 (180195840000 + //(-1346993049600 
+14183769600tt 2 ) + 171953971200;/ 2 + e\ (-1212558336000 + v (1821503692800 
-15274828800tt 2 ) - 124084224000;/ 2 ) + e\ (236763648000 + v (-152134963200 
+1091059200tt 2 ) - 53100748800^) + e 6 t (-9123840000 - 486604800;/+ 5231001600;/ 2 )) 
+ (e t cos u) 7 (-44934912000 + v (424780646400 - 4637001600tt 2 ) - 52218777600;/ 2 
+e? (338950656000 + v (-520930713600 + 4637001600tt 2 ) + 28435968000^) 

+ej (—25774848000— 11146291200;/ + 23782809600;/ 2 )) + (e t cos u) 8 (4561920000 
+// (-50140569600 + 545529600tt 2 ) + 6325862400;/ 2 + e 2 (-47443968000 
+v (65225318400 - 545529600tt 2 ) - 3710361600;/ 2 ) + ej (4561920000 + 243302400;/ 
-2615500800i/ 2 )) + i (e t sin u) (-328593346560 + 26033356800 log (1 - e t cosw) 

-26033356800 log (jpj + 71427240960;/- 11877196800;/ 2 + 2718080000;/ 3 


+e 2 ^962564305920 — 78100070400 log (1 - e t cosw) + 78100070400 log (j- 
+v (-85968921600- 3273177600tt 2 ) + 106786867200;/ 2 - 6687820800;/ 3 ) 
+ef ^—1161059174400 + 78100070400 log (1 - e t cos u) - 78100070400 log ^ 
+v (455668669440 + 4909766400tt 2 ) - 231924172800;/ 2 + 11614080000;/ 3 ) 
+e G t ^483224340480 - 26033356800 log (1 - e t cos u) + 26033356800 log 
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Xo 


+u (-538960911360- 3273177600TT 2 ) + 319443148800;/ 2 - 12469836800;/ 3 ) 

+e® (—22426122240 + 229909708800;/ — 173743948800;/ 2 + 7364044800;/ 3 ) 

+e t 10 (294220800— 9305395200;/ + 23887411200;/ 2 — 1261209600;/ 3 ) + (e t cosu) x 

x ^925176069120— 78100070400 log (1 - e t cosw) + 78100070400 log 

+v (-458849817600+ 8182944000tt 2 ) - 11032934400;/ 2 - 9615590400;/ 3 


+e 2 t ^—2356066314240 + 234300211200 log(l - e t cos u) - 234300211200 log ( 
-300753884160i/- 62224051200^ + 13476403200I/ 3 ) + e 4 t (3159065625600 


x 
V^o 


-234300211200 log (1 - e t cos u) + 234300211200 log 


fe) 


i/(-338253788160 


-81829440007T 2 ) - 32412595200i/ 2 - 11630284800+ 3 ) + e\ (-1380322575360 

+78100070400 log (1 - e t cosm) - 78100070400 log + v (811123860480 

+9819532800tt 2 ) - 319331097600j/ 2 + 5767065600+ 3 ) + e t 8 (48416071680 
-459194112000z/+ 253455436800i/ 2 - 6922828800+ 3 ) + e t 10 (-294220800 
+9305395200*/- 23887411200+*+ 1261209600i/ 3 )) + (e t cos u ) 2 (-901804861440 

+78100070400 log (1 - e t cos u) - 78100070400 log (j^j + v (1567758735360 

-343683648007T 2 ) - 4901222400*/ 2 + 16293504000*/ 3 + e 2 t (1178660367360 

-234300211200 log (1 - e t cos u) + 234300211200 log + v (638531481600 
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-78100070400 log (1 - e t cos u) + 78100070400 log + */ (-12701306880 

-98195328007T 2 ) - 123584793600*/ 2 + 10020940800;/ 3 ) + ef (-31024880640 
+275186073600*/ — 105116083200*/ 2 - 1940428800*/ 3 )) + (e t cos u ) 3 (452697154560 

-26033356800 log (1 - e t cos u) + 26033356800 log + v (-2649367956480 

+58917196800tt 2 ) + 75336499200i/ 2 - 10847257600*/ 3 + e 2 t (1093220843520 

+78100070400 log (1 - e t cos u) - 78100070400 log +1/(5716638720 

-49097664000tt 2 ) - 361501670400i/ 2 + 20455526400*/ 3 ) + e\ (211920660480 

-78100070400 log (1 - e t cosm) + 78100070400 log + v (602037273600 
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+26033356800 log (1 - e t cos u) - 26033356800 log + v (-367892106240 

+3273177600tt 2 ) + 163837209600*/ 2 - 587878400*/ 3 ) + e\ (5034931200 
-459016704001/ + 25404595200*/ 2 + 1499212800*/ 3 )) + (e t cos u ) 4 (-287513978880 
+*/ (2242649057280 - 52370841600tt 2 ) - 56454220800*/ 2 - 971289600*/ 3 
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+*/ (-933238149120+ 24548832000tt 2 ) + 9044275200*/ 2 + 2341017600*/ 3 
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-5811609600*/ 3 ) + e\ (-90193536000+ v (163341803520 + 1636588800tt 2 ) 


25 


—128144793600^ 2 - 3793766400+ 3 ) + e\ (-12266726400- 36156672000^ 

— 1583001600i' 2 - 399667200^ 3 )) + (e t cosu) 6 (-29502489600 + z? (159620889600 
-4909766400 tt 2 ) - 115200000z/ 2 + 81536000+ 3 + e 2 t (-48882355200 + v (6869283840 
+3273177600?r 2 ) + 2940825600^ - 128025600i/ 3 ) + e\ (11003212800 - 50441702400?? 

+30168576000i/ 2 + 1380019200Z/ 3 ) + (1385856000 + 6721920000i/- 422092800z? 2 

-56192000z/ 3 ))))}. (Ale) 
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